/***************************************************************************
 *   Copyright (C) 2007 by Alexandre Masselot,,,   *
 *   alex@alex-laptop   *
 *                                                                         *
 *   This program is free software; you can redistribute it and/or modify  *
 *   it under the terms of the GNU Library General Public License as       *
 *   published by the Free Software Foundation; either version 2 of the    *
 *   License, or (at your option) any later version.                       *
 *                                                                         *
 *   This program is distributed in the hope that it will be useful,       *
 *   but WITHOUT ANY WARRANTY; without even the implied warranty of        *
 *   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the         *
 *   GNU General Public License for more details.                          *
 *                                                                         *
 *   You should have received a copy of the GNU Library General Public     *
 *   License along with this program; if not, write to the                 *
 *   Free Software Foundation, Inc.,                                       *
 *   59 Temple Place - Suite 330, Boston, MA  02111-1307, USA.             *
 ***************************************************************************/
#include "enzymetrypsin.h"
namespace bol{
unsigned char EnzymeTrypsin::m_K=AminoAcid::aa2reduced('K');
unsigned char EnzymeTrypsin::m_R=AminoAcid::aa2reduced('R');
unsigned char EnzymeTrypsin::m_P=AminoAcid::aa2reduced('P');

EnzymeTrypsin::EnzymeTrypsin()
{
  setName("trypsin");
}


EnzymeTrypsin::~EnzymeTrypsin()
{
}

void EnzymeTrypsin::cleavePositions(const Peptide& src, vector<int>&vpos){
  vpos.clear();
  if(src.length()==0) return;
  for(unsigned char *c=src.sequence()+1; *(c+1) !=AminoAcid::CHAR_CTERM_PEPTIDE; c++){
     if((*c==m_K || *c==m_R) && (*(c +1)!=m_P))
        vpos.push_back(c-src.sequence()+1);
  }
}
}




